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The electromagnetic bremsstrahlung spectrum of the dipole which falls along a 
spiral orbit into the Schwarzschild black hole was found. The characteristic features 
of this electromagnetic spectrum can be used for determination of the black hole 
mass. This new way (if implemented) provides higher accuracy of the determination 
of the black hole mass. Also these features in the spectrum can be used for the 
determination of the some physical processes in the black hole magnetosphere and 
in the accretion disk. It is also shown that the asymptotic behavior of this spectrum 
(at high frequencies) is practically independent from the impact parameter of the 
falling dipole. 


I. INTRODUCTION 

As it is well known (see eg DEI) black holes no ’’hair”. Therefore, all electromagnetic 
fields from multipole moments will disappear if the system of charges approaches to the 
black hole horizon. For the point charge at rest the electric held was considered by Linet [3J. 
It has been shown, in particular, that the held of a point charge approaches to the held of 
the charged black hole (with the same charge) when the charge approaches to the horizon. 
Hence it follows that all helds of all electric and magnetic multipole moments should be 
radiated when the charge (or a system of charges or currents) approaches to the black hole 
horizon. 

For accelerated motion of monopole (single charge) the energy loss is determined mainly 
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by bremsstrahlung. This radiation is a dipolar, since the radiated field components are 
inversely proportional to the c 2 (square of the speed of light). 

The existence of dipole radiation is not obvious for a massive dipole falls into a black hole 
(due to the fact that both of the charges of the dipole are moving and accelerating in the 
same direction, and the signs of the charges are opposite). However, the increasing of the 


space-time curvature leads to the emission of the dipole type of radiation (see section IV). 

Earlier, a similar radiation was investigated in many papers - see for example (MO]; We 
do this a slightly different method, which is much more convenient for numerical calculations 
- see [IT, T2| . 

In the paper ra it was found the bremsstrahlung radiation for a charge passing (with a 
constant velocity) through the wormhole. 

In this paper we will solve our problem in general case of non-radial dipole which free 
falling into a Schwarzschild black hole. 


II. LAW OF MOTION OF A FREE-FALLING PARTICLE 


Schwarzschild metric for a nonrotating and uncharged black hole has the form: 

ds 2 = (i_ r ^ c 2 dt 2 dr 2 ^r 2 {d0 2 +sin 2 0 dip 2 ). (1) 

Here: r g = 2 GM/c 2 - the radius of the black hole horizon, G - gravitational constant, M - 
the mass of the black hole. 

Let us remind the law of motion for the test particles in the Schwarzschild gravitational 
field. We suppose that the particle moves in the equatorial plane (9 = 7r/2). Let us write 
the geodesic equation for the particle (see. [HI H5] , §87): 


dui 1 dg ki k j 
—— = — u u 

ds 2 dx l 

Hence for the metric ([!]) and i — 0 (corresponding to the time coordinate ct ) we have the 
integral of motion 1 : uq := e = const (the specific energy the particles); and for i, which 
corresponds coordinates ip we have the integral of motion: = L/(mc) := he = const (here 

L - angular momentum, m - mass and h - the impact parameter of the particle). 


1 In our work, the Latin indices run a series of spatial coordinates and time coordinates, and Greek indices 
run a series of spatial coordinates. 
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We assume that the fall of the particles happens from the initial radius ro and the result 
of this fall is the capture of the particles by a black hole. We also assume that the mass of 
the particle is much smaller than the mass of the black hole: m « M. From equations (J2]) 
and taking into account the identity UiU 1 = 1, we have: 

^ = u° = e/(l -r g /r), (3) 

^ = u r = -Je 2 - (1 - r g /r)( 1 + h 2 e 2 /r 2 ), (4) 

d 4- = u? = -he/r 2 . (5) 

as 

Here r(t) - the current radial coordinate of the particle. 

To find the highest possible value 2 of the impact parameter of the particle h max (in which 
it will be captured by a black hole), we proceed similarly the work [16]. 

The value of h max is defined as the root of the equation [u r ] 2 hr ^ = 0, where this root at 
the same time should be a point of minimum for function [u r ] 2 r y Thus the point of minimum 
r m determined by solution of the equation: 


(« 2 - l)r 2 , + (2 - 1.5£ 2 )r m r g - r 2 = 0. 


( 6 ) 


From the two roots of this equation we need to choose a smaller (which corresponding to the 
plus sign and the minimum distance to the black hole). Then the value of h max is determined 
by the expression: 


7 2 = (e 2 ' r m - r m + r g)r 2 - l 

max 


e 2 (r m -r g ) 


( 7 ) 


At 6 = 1 (particle is at rest at infinity) we have: 


h ( e 1) — 9r* 
u max ^'9' 


( 8 ) 


In the limit e —> oo (for photons), we have: 

} x (e->°o) __ Syfag ( , g'\ 

The minimum valid value for e is determined by the inequality: e 2 > 8/9. 


2 This corresponds to the minimum-possible value of the impact parameter of the particle, in which it still 
will not be captured by a black hole. 
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III. RADIATION OF A CHARGE 


For the covariant 4-vector potential A (for the electromagnetic field), we have two invari¬ 
ants: inv i = AiU 1 and inv 2 = A % A l . In accordance with the gauge invariance of the 4-vector 
potential we choose the calibration A* in the comoving (free-falling) reference system so 
that 3 A 1 := 0. Thus, in the comoving reference frame we have: Aq = q/R, where q - the 
particle charge, R = \r Q b — r\,r 0 b and r- the radius-vectors of the observer and of the charge 
(respectively). Hence, according to the Lorentz transformations, we obtain for the covariant 
spatial components of the potential in the Schwarzschild reference frame near the particle: 


A Q A Q 

Aq = ^u 0 , A 7 = —u 


7‘ 


( 10 ) 


To get the covariant components of the electric field F 0 


gamma 


= d 0 A 7 — 8 7 Aq it is necessary 


to differentiate the expressions (10) with respect to the time coordinate ct 0 b (at the obser¬ 
vation point) and with respect to the spatial coordinates of the radius vector r 0 b (at the 
observation point). Because uq and u v are constants (see section 0 and uq = 0, we have: 


F()r j-fdoZLr. 

lx 


(ii) 


In the expression ( JTT| we have discarded a member with the asymptotic oc 1 /K 2 and al¬ 
lowed only a member with the required asymptotic behavior: oc 1 /R in the field of the 
electromagnetic wave 4 (emw): 


07 


V F07 qc _ 

emw ’ 


pemw _ 

± ± emw ^ ' 


( 12 ) 


The asymptotic behavior oc 1 /R in the electromagnetic field corresponds to the field of 
the wave propagating along the vector at infinity R. In this case, the field asymptotic 
behavior oc 1 /R 2 does not correspond to electromagnetic waves, becouse it does not satisfy 
the conservation of the energy flux through a sphere with the radius r » r g . At the 
same time, at distances which is much larger than the size of the radiating system (ie 
maximum radial coordinates of the particle), wave approximation works and we can talk 


3 In general, the gauge transformation A —► A + dif(x k ) can not be made to vanish by all three spatial 
components Ay, but in the comoving (for incident particle) reference system it is possible to do - see for 
example [13], §65. 

4 Therefore, to find the radiated field components we need to differentiate only Ui. 
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about radiation of photons by the particle. We are interested not only the processes at 
large distances in the wave zone, but (mostly) the processes near the black hole and near 
the particle. When you change locations of the charge (or a system of the charges) in 
Schwarzschild coordinates the corresponding change of electromagnetic field propagates at 
the speed of light. We are interested in the change in electromagnetic field, which corresponds 
to the asymptotic behavior of the wave field (oc 1 /R) and then we will talk about such a 
field - as the radiation field of of the electromagnetic wave with wave vector k l . At the same 
time, we remember that we are considering processes occuring not only in the wave zone. 

To find the Fourier transform of the radiation field we proceed similarly to 5 our work HI]: 


F {w) = - 

r 0 ~/ 9 


[sin (wt oh ) + cos (wtob)] Fo-y dt„b. 


(13) 


Here, the value w = ck$ has the physical meaning of covariance, time component of the wave 
4-vector 6 , A) 0 ~ integral of the motion for the photon (analogue w 0 from (J2])). 

Because the observer with the radial coordinate r ob and 4-velocity U l (relative to predeter¬ 
mined reference system) registers the frequency of the electromagnetic field w := ckiU l , then 
in general the frequency of which registers the fixed observer is equal to w = w f 1 — r g /r ob . 


From the point of view of the stationary observer (at radius r ob ) the formula (13) can be 
rewritten as: 


F {w) ^ J_ 

01 V2 


cos wt ob 


r ob 4 


F 0l dt 0 b‘ 


(14) 


At r ob » r q we have w ~ w. 


In formulas (13) and (14) the point of time t ob is the time of arrival of the wave to the 
observer. At the same time t ob should be carried differentiation (to obtain the field Fo'y from 


5 Here we take into account that up to time t 0 b = 0 the particle at rest (for a stationary observer), therefore 
the emitted fields are absent. Since any function can be represented as a superposition of symmetric and 
antisymmetric parts, then we extrapolating the function F tl to the negative half-time even and odd way 
- in the half-sum we get zero, and each of the parts (even and odd) we use for the Fourier transform. 

6 The values of the contravariant components of the wave 4-vector k l (corresponding to the electromagnetic 
wave which propagating in the plane of the equator) are obtained from expressions (3|5 1 by the replacement 
u l —y eh 1 in the limit e —> oo. 
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potential Ay). Therefore, in expression (14) we can use integration by parts 7 : 


r 0-y ~ 


) _ A^ ob 0) wy/1 - r g /r ob 


+ 


V2 


sm | wt ob J 1 — - j ) A^ ob) dt. 


ob • 


(15) 


In order to determine the time t ob we consider two consecutives near events in the 
Schwarzschild reference system: 

1 ) the field measurement by observer at the time of arrival of the particle at the point with 
the radial coordinate ri; 

2 ) the field measurement by observer at the time of arrival of the particle at the point with 
the radial coordinate r 2 . And rq — r 2 := dr > 0. 

The corresponding total changing of the time dt ob consists of two components: 

1) Period of time 8 A tf, which corresponds to the difference in length of paths 5l± and hZ 2 
for distribution (with speed of light) a retarded potential of the radiated field. Moreover, the 
retarded potential extends from the points with radial coordinates rq and r 2 (respectively) 
to the observer (with the radial coordinate r ob ). If you know the impact parameter for 
the photon hf, then from expressions (J3| and (|4]) at e —> oo (for distribution of retarded 
potential) we get by clock a stationary observer at infinity A tf = /u r )dr / c: 

-y/1 -r g /r ob dr 


At f = 


c( 1 — r g /r) 1 — (1 — r 9 /r)/r^/r 2 


(16) 


2 ) the time required to move the particle from the point r\ to the point r 2 (by clock a 
stationary observer) - see (§ and (§). 

Thus, the sum of these two components gives: 

■a/1 - r g /r ob dr -y/l - r g /r ob dr 


dt nh 


_v_ x ' _ _|_ _v_ z' _ 

c(l - r g /r)yjl - (1 - r g /r)h 2 f /r 2 C ( X “ r g/ r )~ (l ~ r g /r)(l/e 2 + h 2 /r 2 ) 


(17) 


IV. THE RADIATION OF THE DIPOLE 


The held of the retarded potentials of the dipole is determined by the standard way: as 
a superposition of the fields from each of the charges of the dipole. The dipole moment do 


In formula (151 it was taken into account that on the horizon for the dipole we have h 7 9 = 0 - see. (23). 


Hereafter the subscript corresponds to the retarded potential of the field. 
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is also defined by the standard way: 


do q ■ SI — const, (18) 

here SI - constant length of the rigid dipole which is defined in the framework of general 
relativity. To determine the SI we assume that in the reference frame which associated 
with the dipole one of its charges emits a gamma quantum, and the other charge through 
the proper time Sr this gamma quantum are absorbed. Then the square of the interval 
ds 2 between these two events gets (in the comoving and Schwarzschild reference systems, 
respectively): 

ds 2 = c 2 Sr — SI 2 = 0, ds 2 = gijSx l Sx j = 0. (19) 

From here to determine the magnitude SI we use another invariant - see (J2]): 

cSr = Sx l Ui = inv (20) 


We consider three possible variants for the dipole orientation in space: along Sr, along SO 
and along Sip. For the dipole orientation along the increment Sx a (in spatial coordinates) 


we introduce the notation: Sx a := SI • /“. From equations (19) and (20) we obtain for these 
three possible orientations of the dipole: 

(! - r g /r) g _ V 1 - r 9 / r fip _ V 1 ~ r a! r 


r = 


e — u 7 


er 


r + hy /1 — r g /r 


( 21 ) 


Then, in the formula (10), which written for the dipole, it will be a corresponding increment 


Suy for 4-velocity under the parallel translation Uy along Sx a in a curved space-time: 


Ay = -| Su 7 = (dally - r y a ll k ) Sx a = (dally - T k la Uk) 


Q 


d 0 f c 


R 


R 


R 


( 22 ) 


Here Y l kl := \g im ( dig mk + d k g m i ~ d m g k i) - Christoffel symbols. 


The Fourier transform of the electric dipole field at infinity we obtain by substituting 


(22) in the formula (15). 


The retarded potentials (10) and (22) in the leading order are proportional to 
1/c (inversely proportional the speed of light), because in the nonrelativistic limit 
(e —>■ 1, r g /r —>■ 0). The spatial components for 4-velocity inversely proportional to the 
speed of light (w 7 oc v/c), and zero (time) component of the particle 4-velocity in the main 
approximation by c does not depend. Therefore the emitted fields of the charge and the 
dipole in the main approximation are corresponded to the emission of the dipole type radi¬ 
ation (which inversely proportional to the square of the speed of light). 
















V. THE CALCULATION OF THE SPECTRAL DENSITY OF THE DIPOLE 


RADIATION 


We pass to dimensionless units: r g — 1 and c— 1. Let us consider the simplest (in terms 
of computations) dipole orientation 9 : Sx a — 5g ■ 81 ■ / . In this case, the orientation of the 
dipole in spherical coordinates stored over time. Let the observer has the coordinates: 


r ob » r g , 0 ob = 7 t/2 , (p ob = 0 (at the equator plane). Then from (22) we obtain: 

I d 0 f e ru r 


— 


5 ! d 


R 



r) l e 2 r 2 


(23) 


Note that the form of the expression (23) ensures the convergence of the integral (15) 


throughout the range of the radial coordinate r and time t ob (despite the divergence of 
expressions for dt ob near the horizon of the black hole). 

For brevity, we will talk about the photons (instead of the above-mentioned fields of 
retarded potential of the electromagnetic wave), thus all characteristics of the photon is 
determined from the black hole (in the wave zone). Accordingly, for a dipole we introduce 
the quantities subscripts ’V\ and for the held of the retarded potential - 

Denoting: £ := 1/r, from (|3]{5]) we obtain for the total deflection angle of the photon: 

Cd. 


Pf = 


hf d£ 


€06 




+ Pd 


(24) 


During the movement the particle (dipole) emits photons in opposite directions, but we are 
interested in only those photons which are reach the observer, ie at position Pf — p 0 b 0 
and £ 0 6 = l/r ob . In this case, the initial coordinates for photons are: r d , p d (r d ) and time 
t d (r d ) (for each of the photons - its value r d ). To calculate the trajectory of these photons, 
we also need to know their impact parameters hf(r d ). Therefore, to calculate the change of 
the impact parameter dhf (between different photons - during the motion of a particle) we 


obtain from (24) the differential relation: 


dpf(h f ,£d) = di d + ^ dh f + d<p d (i d ) = 0, 


did 


(25) 


9 Here we have designated: 5g - Kronecker symbol. 
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E tot - 0.0011 
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hd = 0.5 h max 
ro = 20 r g 
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E to t “ 0.0027 
max|F< w >| “0.68 
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\/ \AAA^ w, i” a25 

v \ 
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E = 1 
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max|F< w >| - 0.58 
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E tot “ 0.0028 

max|F< w >[ “ 0.66 


(ui) 


1 2 


06 


The parameters h d , ro, e, total 


Figure 1: The spectral density of the dipole radiation F, 
radiated energy E to t (in units d^/(r g R 2 )) and the amplitudes IF^' (in units do/R) for all these 
cases displayed on the panels. The left-hand sides of panels correspond w = 0 ; right-hand sides - 
w = c/r g . 


where 


dh 


d Ff 



did 

\A - hjidi 1 - 

id) 

id 



- [ 

di 


J 

iob 

[i - h}e(i - o] 3/2 ’ 

© we 

have: 



hd 


0- 

-(e- 2 + Ke d )(l 

- id) 


(26) 


(27) 


(28) 


At <fd — 0 we have hf — 0 (the initial moment t d — 0). Tims, from the relation (25) we can 


calculate hf. 
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If we know the value of hf for the photon, we can find the time of the photon propagation 


A tf from the current position of the particle to the observer - see (16): 


id 



£ob 




(29) 


At C, 0 b -4 Owe have: Atf —> oo - since the time of propagation for the photon to an infinitely 
for the distant observer tends to infinity. Therefore, the value must be greater than zero. 


VI. CHARACTERISTIC FEATURES OF THE SPECTRUM 


We will explain the characteristic features of the spectrum of the falling into a black hole 
dipole (for low frequencies w < c/r g ). For this purpose, we according to the most general 
considerations: 

1. At the initial moment a distant observer can detect the electric dipole held, as dipole 
is sufficiently far away from the black hole. 

2. At the final moment (in the limit t Q b —>• oo) the dipole disappears below the horizon 
of the black hole, all electromagnetic fields associated with the dipole must be radiated, so 
a distant observer can not register a dipole held. 

The law of decrease of the held for a distant observer was calculated in the previous 
section. In general, this law can be approximated by any decaying function. In the simplest 


case, we can use the step function: F 0l = F step • ©(£i — £ 0 &), according to the formula (13) 
we have: 


rp( w ) _ p 

-C07 — ^ step ' 


1 — cos(utfi) + sin(rofi 
2 w 


(30) 


If we use for this exponent: Fq 7 = F exp ■ exp(—atf 0 &) ■ 0(£ o & — ti), a > 0, we have: 


F^ ] = F exp ■ exp(-atf 1 ) 


a sin(wtf 1 ) + w cos(utf 1 ) + a cos(wii) — w sin(utf 1 ) 


(31) 


2(ct 2 + w 2 ) 

From the last two formulas it becomes clear why the oscillations in the spectra on Fig. [l] 
and their asymptotic behavior: oc 1 jw. 
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Figure 2: The characteristic shape of the spectrum asymptotics at re —* oo (this expression is 
(35), which squared and multiplied by w 2 ). Panels A and B correspond to the different possible 


.— QQ' ^7’( n ) /q\ £n 

relationships between the coefficients Fq q (0) and y2 0a n\ —~ i* 1 (35) 

n =0 


VII. ASYMPTOTICS OF THE SPECTRUM OF THE RADIATED FIELD AT 


uj —> oo 


To calculate the required asymptotics (at w —> oo - in dimensionless units) we proceed 


similarly to our work Pt The integral in ( Jl5| ) splits into two parts: the first part - from zero 
to the observer’s time t\ and the second - from t\ to infinity (which corresponds to achieve¬ 
ment by the particle the horizon). Moreover, the time t\ is chosen so that the corresponding 
radial coordinate of the particle rq differs from r g by a small amount := rq — r g << r g . 


The potential A 7 in the second part of the integral, according to (23), is equal to 


(d 0 /R) ^j8]_/r g . In addition, according to (17), the element of time dt Q b is equal to 
r g dr/8\. Then in the main approximation by small value di the second part of the 


integral in (15) will be small compared to the first part (oc y/5i). 


Since the main contribution to the integrals (13), (14) and (15) gives the first part (from 


0 to ti), then the second part we can neglect and consider the first integral in (13) in more 
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detail. We decompose the integrand function F 0a (t o b), in a Taylor series with respect to 
dimensionless time t: 


Foa(t) — ^2 


F&h 0 ) • r 


n =0 


n\ 


(32) 


Use the recurrence relation: 


t n nf 

t n [sin(tct) + cos(u>f)] dt — — [sin(wi) — cos(iui)]-/ t n ~ 1 [sin(u>t) — cos(tct)] dt 


t n , . , nt n 1 r . , . n(n — 1) , 2 


= — [sin(wi) — cos(iui)] 


[sin(iui) + cos (tut)] + 


w " w 2 w 2 J 

And we use its mathematical limit (the leading order) at w —)• oo: 
^1 


t n 2 [sin(iui) + cos(tct)] dt 


t n [sin(ret) + cos (wt)] dt 


t™ [sin(ruii) — cos(wii) + 


w 


(33) 


Using this expression in the decomposition (32), for (13) we obtain the asymptotic behavior 
of the spectrum of the radiated field: 


lirn F^ w) « sin ( wt i) -cos(wti) ^^ F^(0)t^ + F 0 q (0 ) 


2w 


<n =0 


n\ 


2 w 


(34) 


This shows that the spectral energy density of the attractor S w of the radiated field (which 
proportional to the square of this expression) is inversely proportional to the square of the 
frequency: S w oc 1/w 2 . 


The expression (34) can be rewritten as: 


lim « ^0) _ A (F ^)g | cos u + r,. 

°“ 2w 2 w\^ n / V 4/ 

\n=0 


(35) 


At w —* oo the shape of this spectrum is almost independent on magnitude of the impact 
parameter h d - see Fig. [2} 

VIII. THE DEPENDENCE OF THE POYNTING VECTOR ON TIME 


Poynting vector P is given by 


10. 


F 2 

p , = r oe 


47rr 2 ’ 


Foe — 


dA.fi 


dt, 


(36) 


ob 


10 Strictly speaking the value of P is not the Poynting vector, but only one of its components (at infinity). 
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Figure 3: Dependencies of the Poynting vectors P := F^g/ [A-nr 2 ) on t ob . The maxima of the curves 
are placed in the center of the panel, the width of the panels - 60r 9 /c, the distance between 
the vertical lines is equal to 15 r g /c, therefore the width of the ’’bells” at half height ~ 10r g /c. 
The amplitude maximum (at different h c j) indicated in units dg /(r' 2 R 2 ) (at r g = 1 and c= 1). 
For each panel, the red curve (with highest amplitude) corresponds to hd = 0, the green curve - 
hd = 0.5h ma x, and the blue curve (with more than one ’’bell”) - hd = 0.99h ma x- 


According to ©> we replace the differentiation with respect dt a b to the differentiation with 
respect r and obtain: 


| — 


doyJ( 1 - r g /r) [r 3 - (r - r g )h)\ [r g r 3 - 2 r g (r - r g )(r 2 /e 2 + h 2 d ) + 2(r - r g ) 2 h 2 d ] 


2 Rr 5 y/l — r g /r 0 b y^r 3 — (r — r g )h 2 + s/r 3 — (r — r g ) (r 2 /e 2 + h d ) 


.(37) 


As it was shown in P)2], and as it is shown in Fig. [3j prohle shape dependent on the Pointing 
vector on time t 0 b is practically independent from the dipole orbital parameters (the initial 
radius of falling). And the prohle width at half maximum A t 0 b directly related with the mass 
of the black hole: A t a b ~ 10r s /c. In general case, non-radial dipole falling into Schwarzschild 
black hole the prohle form P(t ob ) becomes dependence on e and r 0 at hd > 0 (see Fig. [3]). 
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At hd ~ 1 in the function P(t 0 b) appears several local maxima, although the maximum 
amplitude is reduced and at the same time slightly reduced the total radiated energy - see 

Fig-13 


IX. DISCUSSION AND CONCLUSIONS 


It is seen from the results of the numerical calculations (see. Fig. Q, the main features of 
the spectrum are also saved for hd > 0 (asymptotic behavior and an increase in the frequency 
of the oscillations with increasing ro). In this case, all the spectral dependence (on the value 
of hd) is detectable only at small w <~ 1, and at the asymptotic behavior (w —> oo) (as 


shown in Section VII). The behavior of the spectra does not vary. 


The dependence of the amplitude and the total radiated energy on the value hd is weak, 
therefore almost all of the results and conclusions obtained by us for the case hd = 0, are 
also applicable to the general case hd > 0. We list here all these results (see P2|): 

1. In the case of supervisory detection of such spectra (or rather their asymptotics) we 
have a real possibility a new and independent way to determine the main characteristics 
of a black hole - its mass. In addition, probably also we will be able to determine some 
properties of the magnetized accreating matter. 

2. The main problem, would be too weak energy flow (for measuring on the available 
radio frequencies) for observation of such process (dipole falling into a black hole). 

3. In addition, the falling of the strongly magnetized compact bodies onto black holes 
are apparently very rare events, which also represents a serious obstacle for observations. 

4. It is also likely that the scan of the entire range of the radio waves on the subject the 
search of such spectra lead to the discovery of new black holes in our galaxy and studying 
of their properties. 
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